Striped antiferromagnetism and electronic structures of SrFeAsF 
and their implications 
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Abstract. - We investigate structural, magnetic, and electronic properties of SrFeAsF as a new 
parent for superconductors using a state-of-the-art density-functional theory method. Calculated 
results show that the striped antiferromagnetic order is the magnetic ground state in the Fe 
layer and the interlayer magnetic interaction is tiny. Calculated As and Sr positions are in good 
agreement with experiment. There are only two quasi- two-dimensional bands near the Fermi 
level. The valence charge is mainly in the Fe and F layers, and the magnetic moment is confined 
to the Fe atoms. All the spin couplings within the Fe layer are antiferromagnetic due to the 
superexchange through the nearest As atoms. These results, with the record-equaling phase- 
transition temperature in the latest Sm-doped SrFeAsF, show that the SrFeAsF, sharing the same 
structure with LaFeAsO, is promising for achieving better superconductors. 



Introduction. — Fe-based superconductors attract 
more and more attention since superconductivity was 
found in doped LaFePO [1]. The advent of superconduct- 
ing F-doped LaFeAsO stimulates a world-wide campaign 
for more and better Fe-based superconductors [2] . Replac- 
ing La by other lanthanides or doping with F has yielded 
more superconductors, and higher phase-transition tem- 
peratures (Tc) have been achieved in some of them [3-5] . 
Furthermore, much more superconducting materials were 
found by using other dopants and other parent com- 
pounds [6-9] . Even a FeSe can be made superconducting 
by applying high pressure [10,11]. Now there are three 
series of FeAs-based superconductors: i?FeAsO {R: lan- 
thanide elements), AFe2As2 {A: alkaline-earth elements), 
and LiFeAs. So far, the highest Tc is 55-56 K in the case 
of doped SmFeAsO [5]. Their structural, magnetic, elec- 
tronic properties are intensively investigated and the mi- 
croscopic mechanism for the superconductivity in these 
materials has been explored [12-23]. Very recently, su- 
perconductivity was found in Co and La doped SrFeAsF 
materials [24-27]. SrFeAsF has the same crystal structure 
as i?FeAsO and similar magnetic instability, but does not 
include any lanthanide [28-30]. It has been found that 
Sm-doped SrFeAsF can become superconducting at 56 K, 



and higher transition temperatures should be reached with 
appropriate dopants [25-27]. Because some magnetic fluc- 
tuations are believed to mediate the superconductivity in 
FaAs-based materials, it is highly desirable to investigate 
the magnetic orders, electronic structures, and magnetic 
properties of the parent compound SrFeAsF. 

Here, we use a state-of-the-art density-functional theory 
(DFT) method to investigate the structural, electronic, 
and magnetic properties of the SrFeAsF. Our total energy 
results show that the striped antiferromagnetic (AF) or- 
der, the same as that of LaFeAsO, is the magnetic ground 
state in the Fe layer, and the interlayer magnetic inter- 
action is tiny. Our calculated position parameters of As 
and Sr are in good agreement with experiment. The elec- 
tronic band result shows that there are only two quasi-two- 
dimensional (quasi-2D) bands near the Fermi level. Our 
charge and magnetization density analysis shows that the 
valence charge is mainly distributed in the Fe and F layers, 
and the magnetic moment is confined to the Fe layer. The 
spin couplings within the Fe layer are AF due to superex- 
change through the nearest As atoms. The real moment 
should be much smaller than the DFT value because of 
quantum spin fluctuations. More detailed results will be 
presented in the following. 



p-1 



Li-Fang Zhu and Bang-Gui Liu 






u i 








"s 








(c) 



Fig. 1: (color online). The checkerboard structure [AFl, (a) 
the stripe one [AF2, (b)], and the first Brillouin zone of the 
AF2 structure (c). The largest ball (cyan or white) denotes 
Sr, the medium (yellow or white) As, the smallest (magenta or 
gray) F, and the ball with arrow Fe (red or gray). The arrow 
denotes the spin orientation. 



Computational detaiL — SrFeAsF usually assumes 
the tetragonal AsCrSiZr (tP8) structure with space group 
P4/nnim (No. 129) at high temperature. It transits to an 
orthorhombic distorted phase at 185 K, and furthermore, 
there is a magnetic instability at 175 K [28,29]. A similar 
behavior has been observed in LaFeAsO [2, 12, 13]. We 
shall investigate by self-consistent calculations the possi- 
ble magnetic structures of the two SrFeAsF phases shown 
in fig. [T] Figure [1^ shows the structure with the famil- 
iar checkerboard AF (AFl) order and fig. [T]d that with a 
striped AF (AF2) order. The interlayer magnetic interac- 
tion can be ferromagnetic (FM) or AF. 

All the calculations are done using a full-potential 
linearized-augmented-plane-wave (FLAPW) method 
within the density-functional theory [31, 32], as im- 
plemented in the package WIEN2K [33, 34]. The 
generalized gradient approximation (GGA) to the 
exchange-correlation potential is used for the presented 
results [35], and local-density-approximation (LDA) 
calculations are done for comparison [36]. Full relativistic 
effects are calculated for core states, and the scalar 
relativistic approximation is used for valence states. The 
spin-orbit coupling [37] is neglected because it has little 
effect on the system. We use 400 k points in the first 
Brillouin zone for the two AF structures. We make the 
harmonic expansion up to /max— 10 inside the atomic 
spheres. The radii of the muffin-tin spheres of Sr, Fe, 
As, and F are 2.3, 2.1, 2.2 and 2.0 atomic unit (a.u.), 
respectively. i?mtX^'^max is set to 8.0. The self-consistent 
calculations are considered to be converged only when 
the integrated charge difference per formula unit between 
input and output charge density is less than 0.0001. 

Results and discussion. — We use the experimental 
lattice constants (a,c) — (3.9930A,8.9546A) for the tetrag- 
onal structure, and (a,6,c) = (5.6155A,5.6602A,8.9173A) 
for the orthorhombic structure [28,29]. As for the internal 
position parameters of As and Sr atoms, we use (zAsj-^Sr) 
= (0.6527,0.1598) and (0.6494,0.1635) as input and opti- 
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Fig. 2: (color online). Total density of states (black thick solid, 
in states/eV per formula unit) and those projected in atomic 
spheres of F (pink or gray dash), Fel (blue or gray dot), Fe2 
(blue or gray short dash), and As (magenta or gray dash dot) 
and the interstitial region (black thin sold) of the AF2 struc- 
ture. 



mize them in terms of forces standard (2 mRy/a.u.) for 
both of the structures [28,29]. We obtain (zAs,^Sr) = 
(0.6444,0.1634) and (0.6475,0.1637) for the AFl and AF2 
structures, respectively. For the AF2 structure, our GGA 
result, (0.6475,0.1637), is in good agreement with experi- 
ment [28,29], in contrast to that from non-magnetic calcu- 
lations [38,39] and (0.6357,0.1624) of an LDA calculation 
of ours. The magnetic moment in the Fe sphere is 1.65 
and 1.97^B for the AFl and AF2 phases. The total Fe 
moment in the AF2 phase is approximately 2fj,B, twice as 
large as an LDA result of ours. The parameters zas and zsr 
and the magnetization density distributions do not change 
when we switch the interlayer magnetic interaction from 
FM to AF. The AF2 phase is lower by 80 meV in total 
energy per formula unit than the AFl phase. This means 
that the striped AF2 phase is the ground-state phase, in 
agreement with experiment [28,29]. The Fe spins align 
parallel along the a (shorter) axis and antiparallel along 
the b axis, which is the same as LaFeAsO [6,17,40]. Sr- 
FeAsF has the striped AF order as its ground state, the 
same as LaFeAsO [16,17,22,23]. The energy difference 
of 80 meV of the AFl phase is also comparable with the 
75 meV for the checkerboard AF phase of LaFeAsO [16]. 
The moment of the magnetic ground state phase is ap- 
proximately the same as that of LaFeAsO [16,22,23]. 

The total spin-dependent density of states (DOS) and 
those projected in the muffin-tin spheres of F, Fel, Fe2, 
and As atom and in the interstitial region of the AF2 
structure are presented in fig. O The filled states between 
-6.9 and eV are the valence states from Fe-3d^4s^, As- 
4p'^, Sr-5s^, and F-2p^ orbitals. The semi-core states such 
as As-4s^ are lower than -10.7 eV. There is no energy gap 
in the energy window presented, but two pseudo-gaps are 
visible, at the Fermi level and -3.85 eV. The DOS be- 
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Fig. 3: Spin-dependent band structure of the AF2 structure 
along representative high-symmetry lines in the first Brillouin 
zone. The left part is that with Fel spin-up character shown 
and the right part with Fel spin-down character shown. The 
band consists of dots, where the bigger the dot, the stronger 
Fel character is at that point. 



tween -6.9 and -5.2 eV, having symmetry between the two 
spin channels, mainly comes from F-2p states, and the 
Sr states almost disappear from the energy window pre- 
sented. The As DOS is also almost symmetrical. Concern- 
ing spin, it is distributed between -5.2 and 2.5 eV, but its 
largest part is between -3.8 and -2.2 eV. It is clear that 
the Fe states are spin-split. Actually, we have two differ- 
ent Fe atoms: Fel and Fe2. Their spins are antiparallel. 
The spin-down part of the Fel DOS is mainly distributed 
between -1 and 2.5 eV and the spin- up part between -3.8 
and -1 eV. The Fe2 DOS is equivalent to the Fel DOS 
with the two spin channels interchanged. One feature is 
a substantial DOS contribution from the interstitial re- 
gion, which reflects substantial covalence between Fe and 
As. At the Fermi level, the DOS calculated without spin 
polarization is substantially larger than that of our spin- 
polarized calculation [38,39]. Our DOS results are similar 
to spin-polarized GGA results of LaFeAsO [16]. 

In fig. [3] we show the spin-dependent band structure 
with Fel character. The bands have very little dispersion 
along the z axis. It is surprising that the band struc- 
ture near the Fermi level is very simple, consisting of only 
two bands at the Fermi level. The band that has the 
maximum along F-Z originates from the Fel A^y state, 
and the other the Fel d^2_y2 one. These two bands near 
the Fermi level have a much larger weight in the spin- 
down than in the spin-up channel. This implies that the 
electronic structure near the Fermi level is quasi-2D and 
spin-split. This special property is reasonable because the 
Sr atom is strongly ionic and the Fe moments align an- 
tiferromagneticaly along the y axis. The band structure 
near the Fermi level is similar to that of LaFeAsO [16]. It 
is interesting that the band feature near the Fermi level 
along T-Y is consistent with that obtained with a phe- 
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Fig. 4: (color online). The charge density distributions in the 
two Fe-As-Fe planes with the Fe atoms along the x axis (a) 
and the y axis (e), the two different Sr-As-Sr planes (c) and 
(g), and the two F-Sr-F planes with the F atoms along the 
X axis (d) and the y axis (h); and the magnetization density 
distributions in the same two Fe-As-Fe planes (b) and (f ) . The 
charge density increment is 0.003e/a.u.^ and the magnetization 
one ±0.003/^s/a.u.^. The red or gray dot in the charge plots 
labels the smallest density. In the magnetization plots, the red 
or gray letters 'z', 'p', and 'n' mean that the magnetization 
density therein is zero, positive, and negative, respectively. 



nomenological two-band model [41]. Actually, the detailed 
band structure near the Fermi level depends sensitively on 
the structural parameters (^Asi^Sr) and the magnetic mo- 
ment [16,22,23]. When zas becomes smaller, the moment 
is correspondingly smaller and the band structure near the 
Fermi level looks more like that of non-magnetic calcula- 
tions [18,22,38,39]. 

In fig. m we present the charge and magnetization den- 
sity distributions in the planes defined by the two sym- 
metrical bonds with the same vertex atom. The panels 
(a) and (e) show the charge density distributions in the 
two planes defined by the two different Fe-As-Fe chains, 
respectively, and the (b) and (f ) the corresponding magne- 
tization density distributions. One of the planes includes 
the X axis and the other the y axis. The panels (c) and 
(g) show the charge density distributions in the two planes 
defined by the two different Sr-As-Sr chains. The panels 

(d) and (h) show the charge density distributions in the 
two planes defined by the two different F-Sr-F chains, re- 
spectively. The red or gray dot in the charge plots labels 
the smallest charge density in e/a.u.^ (e is the electron 
charge): 0.007 in (a), 0.003 in (c), 0.003 in (d), 0.006 in 

(e) , 0.003 in (g), and 0.003 in (h). The lowest contours rep- 
resent the charge density values in e/a.u.'^: 0.004, 0.007, 
0.010, and so on. In the magnetization plots in (b) and 

(f) , the red or gray letters 'z', 'p', and 'n' mean that the 
magnetization density in the region is zero, positive, and 
negative, respectively. The magnetization density values 
(in /^s/a.u.^) of the contours near the zero are ±0.003, 
±0.006, and so on. The magnetization density vanishes in 
these planes at the same places as in the panels (c), (d), 

(g) , and (h), and thus this is not shown again. 

In order to investigate the charge and magnetization 
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Fig. 5: (color online). The charge density distributions in the 
Fe plane (a), the Sr plane (c), the As plane (d), and the F 
plane (f); and the magnetization density distributions in the 
Fe plane (b) and the As plane (e). The red or gray dot labels 
the smallest density. The red or gray letters 'z', 'p', and 'n' 
and the density increments are the same as in Fig. 4. 

density distributions as a whole, we present them by 
atomic layer in fig. [51 The panels (a) and (d) show the 
charge density distributions in the Fe and As planes, and 
the panels (b) and (e) the magnetization density distribu- 
tions in the same planes. The panels (c) and (f) show the 
charge density distributions in the Sr and F planes, and 
the magnetization density in these two planes is near zero 
and thus not presented. The red or gray dots denote the 
smallest charge density in e/a.u.'^: 0.024 in (a), 0.0025 in 
(c), 0.006 in (d), and 0.0025 in (f). The letters 'z', 'p', 
and 'n' imply the same as in fig. [H It is clear that the 
spin is mainly confined in the Fe layer and the magnetiza- 
tion density is nearly zero in the F layer and the adjacent 
two Sr layers. Therefore, the interlayer spin interaction 
must be tiny, which in return supports our computational 
model. 

In addition, we have done the same DFT calculations 
assuming an AF order along the z axis. The calculated 
total energies and position parameters of As and Sr are 
the same as those presented above. This is caused by the 
strong ionicity of the Sr layers and the zero magnetiza- 
tion density in the Sr and F layers. Actually, one cannot 
determine the magnetic order along the z axis by using 
DFT calculations because the interlayer magnetic interac- 
tion is too weak. As for the spin order in the Fe layer, 
the spins align antiferromagnetically along the y axis but 
ferromagnetically along the x axis. The nearest spin cou- 
pling constants along the x and y axes, Jx and Jj,, are 
AF. There is an AF couphng J' (> Ja:/2) between the 
next nearest Fe spins. Other inter-spin couplings are sub- 
stantially smaller [42] . These three coupling constants are 
determined by the super-exchange through the bridging 
As atoms and other factors [40]. The small difference 



S = Jy — Jx is caused by the structural distortion. As 
for the moment per Fe atom, our GGA and LDA results 
are nearly 2^b and much larger than the experimen- 
tal value ~ 0.3/iB [30]. This situation is the same as that 
of LaFeAsO whose experimental moment per Fe atom is 
0.25 ^ 0.36 [12,13], much smaller than 1.5 ~ 2.3 hb 
from GGA calculations [22,23]. For LaFeAsO, the large 
discrepancies can be reduced by using LDA rather than 
GGA in optimizing the crystal structure and calculating 
the moment, but this scheme leads to a large discrepancy 
of the calculated As position parameter from the exper- 
imental result [17,22,23]. On the other hand, spin-orbit 
coupling, monoclinic distortion, and p-d hybridization are 
used to obtain smaller value for the moment [21]. In fact, 
this situation, however, has not been settled even in the 
case of LaFeAsO. We believe that quantum many-body 
effects may play some important roles in determining the 
actual moment. For SrFeAsF, the actual magnetic mo- 
ment is also much smaller than the GGA value and we also 
attribute this discrepancy to quantum many-body effects 
of the spin fluctuations. More investigations are in need 
to solve this moment problem without leading to other 
discrepancies. 

Conclusion. — In summary, we have investigated the 
structural, magnetic, and electronic properties of SrFeAsF 
using the state-of-the-art DFT method. Our total energy 
results show that the striped AF order is the magnetic 
ground state in the Fe layer and the interlayer magnetic 
interaction is tiny. The calculated As and Sr positions 
are in good agreement with experiment. There are only 
two quasi-2D bands near the Fermi level. The valence 
charge is mainly distributed in the Fe and F layers, and 
the magnetic moment is confined to the Fe atoms. All 
the intra-layer spin couplings are AF due to the superex- 
change through As atoms. SrFeAsF shares the main fea- 
tures of the structural, electronic, and magnetic proper- 
ties with LaFeAsO. The actual moment is smaller than 
our value, which is the same as the situation for LaFeAsO 
and should be solved by considering quantum many-body 
effects. Because the Sm-doped SrFeAsF has equaled the 
Tc record (56 K) of the FeAs-based superconductors, it is 
believed that higher Tc will be realized in the SrFeAsF se- 
ries. These results are useful for understanding the struc- 
tural, electronic, and magnetic properties of SrFeAsF and 
should have implications achieving better superconductors 
by appropriate doping. 
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